General Relativistic Counterpart Research Articles

Tidal Deformability of Neutron Stars in Scalar-tensor Theories of Gravity

Gravitational waves from compact binary coalescences are valuable for testing theories of gravity in the strong field regime. By measuring neutron star tidal deformability using gravitational waves from binary neutron stars, stringent constraints were placed on the equation of state of matter at extreme densities. Tidal Love numbers in alternative theories of gravity may differ significantly from their general relativistic counterparts. Understanding exactly how the tidal Love numbers change will enable scientists to untangle physics beyond general relativity from the uncertainty in the equation of state measurement. In this work, we explicitly calculate the fully relativistic l ≥ 2 tidal Love numbers for neutron stars in scalar-tensor theories of gravitation. We use several realistic equations of state to explore how the mass, radius, and tidal deformability relations differ from those of general relativity. We find that tidal Love numbers and tidal deformabilities can differ significantly from those in general relativity in certain regimes. The electric tidal deformability can differ by ∼200%, and the magnetic tidal deformability differs by ∼300%. These deviations occur at large compactnesses (C = M/r ≳ 0.2) and vary slightly depending on the equation of state. This difference suggests that using the tidal Love numbers from general relativity could lead to significant errors in tests of general relativity using the gravitational waves from binary neutron star and neutron star black hole mergers.

Shadows and photon rings of regular black holes and geonic horizonless compact objects

Abstract The optical appearance of a body compact enough to feature an unstable bound orbit, when surrounded by an accretion disk, is expected to be dominated by a luminous ring of radiation enclosing a central brightness depression typically known as the shadow. Despite observational limitations, the rough details of this picture have been now confirmed by the results of the Event Horizon Telescope (EHT) Collaboration on the imaging of the M87 and Milky Way supermassive central objects. However, the precise characterization of both features—ring and shadow—depends on the interaction between the background geometry and the accretion disk, thus being a fertile playground to test our theories on the nature of compact objects and the gravitational field itself in the strong-field regime. In this work we use both features in order to test a continuous family of solutions interpolating between regular black holes and horizonless compact objects, which arise within the Eddington-inspired Born–Infeld theory of gravity, a viable extension of Einstein’s general relativity (GR). To this end we consider seven distinctive classes of such configurations (five black holes and two traversable wormholes) and study their optical appearances under illumination by a geometrically and optically thin accretion disk, emitting monochromatically with three analytic intensity profiles previously suggested in the literature. We build such images and consider the sub-ring structure created by light rays crossing the disk more than once and existing on top of the main ring of radiation. We discuss in detail the modifications as compared to their GR counterparts, the Lyapunov exponents of unstable nearly-bound orbits, as well as the differences between black hole and traversable wormholes for the three intensity profiles. In addition we use the claim by the EHT Collaboration on the radius of the bright ring acting (under proper calibrations) as a proxy for the radius of the shadow itself to explore the parameter space of our solutions compatible with such a result.

Revisiting Vainshtein screening for fast N-body simulations

We revisit a method to incorporate the Vainshtein screening mechanism in N-body simulations proposed by R. Scoccimarro in [1]. We further extend this method to cover a subset of Horndeski theories that evade the bound on the speed of gravitational waves set by the binary neutron star merger GW170817. The procedure consists of the computation of an effective gravitational coupling that is time and scale dependent, G eff (k,z), where the scale dependence will incorporate the screening of the fifth-force. This is a fast procedure that when contrasted to the alternative of solving the full equation of motion for the scalar field inside N-body codes, reduces considerably the computational time and complexity required to run simulations. To test the validity of this approach in the non-linear regime, we have implemented it in a COmoving Lagrangian Approximation (COLA) N-body code, and ran simulations for two gravity models that have full N-body simulation outputs available in the literature, nDGP and Cubic Galileon. We validate the combination of the COLA method with this implementation of the Vainshtein mechanism with full N-body simulations for predicting the boost function: the ratio between the modified gravity non-linear matter power spectrum and its General Relativity counterpart. This quantity is of great importance for building emulators in beyond-ΛCDM models, and we find that the method described in this work has an agreement of below 2% for scales down to k ≈ 3h/Mpc with respect to full N-body simulations.

Perturbations of Spinning Black Holes beyond General Relativity: Modified Teukolsky Equation

The detection of gravitational waves from compact binary mergers by the LIGO/Virgo collaboration has, for the first time, allowed for tests of relativistic gravity in the strong, dynamical and nonlinear regime. Outside Einstein's relativity, spinning black holes may be different from their general relativistic counterparts, and their merger may then lead to a modified ringdown. We study the latter and, for the first time, derive a modified Teukolsky equation, i.e., a set of linear, decoupled differential equations that describe dynamical perturbations of non-Kerr black holes for the radiative Newman-Penrose scalars $\Psi_0$ and $\Psi_4$. We first focus on non-Ricci-flat, Petrov type D black hole backgrounds in modified gravity, and derive the modified Teukolsky equation through direct decoupling and through a new approach, proposed by Chandrasekhar, that uses certain gauge conditions. We then extend this analysis to non-Ricci-flat, Petrov type I black hole backgrounds in modified gravity, assuming they can be treated as a linear perturbation of Petrov type D, black hole backgrounds in GR by generalizing Chandrasekhar's approach, and derive the decoupled modified Teukolsky equation. Our work lays the foundation to study the gravitational waves emitted in the ringdown phase of black hole coalescence in modified gravity for black holes of any spin. Our work can also be extended to compute gravitational waves emitted by extreme mass-ratio binary inspirals in modified gravity.

Shadow revisiting and weak gravitational lensing with Chern-Simons modification

Dynamical Chern-Simons (dCS) gravity has been attracting plenty of attentions due to the fact that it is a parity-violating modified theory of gravity that corresponds to a well-posed effective field theory in weak coupling approximation. In particular, a rotating black hole in dCS gravity is in contrast to the general relativistic counterparts. In this paper, we revisit the shadow of analytical rotating black hole spacetime in dCS modified gravity, based on which we study the shadow observables, discuss the constraint on the model parameters from the Event Horizon Telescope (EHT) observations, and analyze the real part of quasi-normal modes (QNMs) in the eikonal limit. In addition, we explore the deflection angle in weak gravitational field limit with the use of Gauss-Bonnet theorem. We find that the shadow related physics and the weak gravitational lensing effect are significantly influenced by the CS coupling, which could provide theoretical predictions for a future test of the dCS theory with EHT observations.

Compact stellar structures in Weyl geometric gravity

We consider the structure and physical properties of specific classes of neutron, quark, and Bose-Einstein condensate stars in the conformally invariant Weyl geometric gravity theory. The basic theory is derived from the simplest conformally invariant action, constructed, in Weyl geometry, from the square of the Weyl scalar, the strength of the Weyl vector, and a matter term, respectively. The action is linearized in the Weyl scalar by introducing an auxiliary scalar field. To keep the theory conformally invariant the trace condition is imposed on the matter energy-momentum tensor. The field equations are derived by varying the action with respect to the metric tensor, Weyl vector field and scalar field. By adopting a static spherically symmetric interior geometry, we obtain the field equations, describing the structure and properties of stellar objects in Weyl geometric gravity. The solutions of the field equations are obtained numerically, for different equations of state of the neutron and quark matter. More specifically, constant density stellar models, and models described by the stiff fluid, radiation fluid, quark bag model, and Bose-Einstein condensate equations of state are explicitly constructed numerically in both general relativity and Weyl geometric gravity, thus allowing an in depth comparison between the predictions of these two gravitational theories. As a general result it turns out that for all the considered equations of state, Weyl geometric gravity stars are more massive than their general relativistic counterparts. As a possible astrophysical application of the obtained results we suggest that the recently observed neutron stars, with masses in the range of $2{M}_{\ensuremath{\bigodot}}$ and $3{M}_{\ensuremath{\bigodot}}$, respectively, could be in fact conformally invariant Weyl geometric neutron or quark stars.

Observational signatures of Schwarzschild-MOG black holes in scalar–tensor–vector gravity: images of the accretion disk

The present paper contains the investigations with respect to the optical appearance, redshift distribution, and observed energy flux of the thin accretion disk described by the Novikov-Thorne model in modified gravity (MOG) proposed by Moffat. The results show that the gravity coupling parameter alpha plays an indispensable role in the size of direct and secondary images of the accretion disk, while the image profile is mainly affected by the observational angles. In addition, an increase of the parameter alpha leads to that of the redshift factor, but a decrease in the radiation energy flux of the accretion disk. These results are due to the enhancement of the parameter alpha that strengthens the black hole’s gravitational effects. By comparing the accretion disk images of the Schwarzschild-MOG black hole, the Schwarzschild black hole, and the Reissner–Nordström black hole, we point out that the Schwarzschild-MOG black hole exhibits the largest inner shadow and the lowest luminosity. This finding provides a potential way to distinguish Schwarzschild-MOG black holes from their counterparts in General Relativity.

Dressed black holes in the new tensor–vector–scalar theory

As incarnations of gravity in its prime, black holes are arguably the best target for us to demystify gravity. Keeping in mind the prominent role black holes play in gravitational wave astronomy, it becomes a must for a theory to possess black hole solutions with only measurable departures from their general relativity counterparts. In this paper, we present black holes in a tensor-vector-scalar representation of relativistic modified Newtonian dynamics. We find that the theory allows Schwarzschild and nearly-Schwarzschild black holes as solutions, while the nontrivial scalar and vector fields generally diverge at the event horizon. Whether this is a physical pathology or not poses a challenge for these solutions, and by extension, the model. However, even if it is, this pathology could be overcome when the black hole hair vanishes.

Neutron stars in degenerate higher-order scalar-tensor theories

We study neutron star configurations in a simple shift-symmetric subfamily of degenerate higher-order scalar-tensor (DHOST) theories, whose deviations from General Relativity (GR) are characterized by a single parameter. We compute the radial profiles of neutron stars in these theories of modified gravity, using several realistic equations of state for the neutron star matter. We find neutron stars with masses and radii significantly larger than their GR counterparts. We then consider slowly-rotating solutions and determine the relation between the dimensionless moment of inertia and the compactness, relation that has the property to be almost insensitive to the equation of state.

Cosmic censorship and charged radiation in second order Lovelock gravity

The conditions for naked singularity formation are considered for a radiating metric of Boulware–Deser type within an electromagnetic field in second order Lovelock (or Einstein–Gauss–Bonnet) gravity. The spacetime metric remains real only up to certain maximum charge contribution. This differs from general relativity. Beyond a certain maximal charge, there exists no real and physical spacetime since the metric becomes complex. We establish that, under certain parameters and for specific values of the mass function and charge contribution, this branch singularity is indeed a naked singularity. This is in contrast to the neutral case where the spacetime metric is always real for a positive mass function, and further, a weak, initially naked singularity always occurs before it becomes covered by an event horizon for all future time. We highlight that both neutral and charged collapse under gravity in Einstein–Gauss–Bonnet gravity differ significantly to their general relativistic counterparts.

Regular black holes without mass inflation instability

Generic models of regular black holes have separate outer and inner horizons, both with nonzero surface gravity. It has been shown that a nonzero inner horizon surface gravity results in exponential instability at the inner horizon controlled by this parameter. This phenomenon takes the name of “mass inflation instability”, and its presence has put in question the physical viability of regular black holes as alternatives to their (singular) general relativity counterparts. In this paper, we show that it is possible to make the inner horizon surface gravity vanish, while maintaining the separation between horizons, and a non-zero outer horizon surface gravity. We construct specific geometries satisfying these requirements, and analyze their behavior under different kinds of perturbations, showing that the exponential growth characteristic of mass inflation instability is not present for these geometries. These “inner-extremal” regular black holes are thereby better behaved than singular black holes and generic regular black holes, thus providing a well-motivated alternative of interest for fundamental and phenomenological studies.

Fingerprints of modified gravity on galaxies in voids

ABSTRACT We search for detectable signatures of f(R) gravity and its chameleon screening mechanism in the baryonic and dark matter (DM) properties of simulated void galaxies. The enhancement of the gravitational acceleration can have a meaningful impact on the scaling relations as well as on the halo morphology. The galaxy rotational velocity field (calculated with the velocity of the gas disc and the acceleration fields) deviates from the typical values of the Tully–Fisher Relation in General Relativity (GR). For a given stellar mass, f(R) gravity tends to produce greater maximum velocities. On the other hand, the mass in haloes in f(R) gravity is more concentrated than their counterparts in GR. This trend changes when the concentration is calculated with the dynamical density profile, which takes into account the unscreened outer regions of the halo. Stellar discs interact with the overall potential well in the central regions, modifying the morphology of the screening regions and reshaping them. We find a trend for galaxies with a more dominant stellar disc to deviate further from round screening regions. We find that small haloes are less triaxial and more round in f(R) than their GR counterparts. The difference between halo morphology becomes smaller in f(R) haloes whose inner regions are screened. These results suggest possible observables that could unveil modified gravity effects on galaxies in voids in future cosmological tests of gravity.

Anisotropic quark stars in f(R) = R 1+ϵ gravity

Within the metric formalism of f(R) theories of gravity, where R is the Ricci scalar, we study the hydrostatic equilibrium structure of compact stars with the inclusion of anisotropic pressure. In particular, we focus on the f(R) = R 1+ϵ model and we examine small deviations from general relativity for |ϵ| ≪ 1. A suitable definition of mass function is explicitly formulated from the field equations and the value of the Ricci scalar at the center of each star is chosen such that it satisfies the asymptotic flatness requirement. We find that both the mass and the radius of a compact star are larger with respect to the general relativistic counterpart. Furthermore, we remark that the substantial changes due to anisotropy occur mainly in the high-central-density region.

Constraining spontaneous black hole scalarization in scalar-tensor-Gauss-Bonnet theories with current gravitational-wave data

We examine the constraining power of current gravitational-wave data on scalar-tensor-Gauss-Bonnet theories that allow for the spontaneous scalarization of black holes. In the fiducial model that we consider, a slowly rotating black hole must scalarize if its size is comparable to the new length scale $\ensuremath{\lambda}$ that the theory introduces, although rapidly rotating black holes of any mass are effectively indistinguishable from their counterparts in general relativity. With this in mind, we use the gravitational-wave event GW190814---whose primary black hole has a spin that is bounded to be small, and whose signal shows no evidence of a scalarized primary---to rule out a narrow region of the parameter space. In particular, we find that values of $\ensuremath{\lambda}\ensuremath{\in}[56,96]\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ are strongly disfavored with a Bayes factor of 0.1 or less. We also include a second event, GW151226, in our analysis to illustrate what information can be extracted when the spins of both components are poorly measured.

Rigid rotation in GR and a generalization of the virial theorem for gravitomagnetism

In this work we study the properties of rigidly rotating neutral dust solutions in general relativity. This class of solutions gained relevance recently due to applications to the dynamics of spiral galaxies. We show that this class could be interpreted as a “rigid body” in general relativity and we analyze the different properties respect to the rigidly rotating disk in special relativity: for example, the general relativistic counterpart shows no Doppler effect for a light signal emitted and received from any two points at rest respect to the “rigid body”. This effect can be important to test the validity of the assumed model for our galaxy. In the second part we approach the problem from a low energy expansion perspective and we write down a generalization of the virial theorem for stationary spacetimes. The non-Newtonian contributions can lead to a re-weighting of dark matter in galaxies.

Black hole sensitivities in Einstein-scalar-Gauss-Bonnet gravity

The post-Newtonian dynamics of black hole binaries in Einstein-scalar-Gauss-Bonnet theories of gravity depends on the so-called "sensitivities", quantities which characterize a black hole's adiabatic response to the time-dependent scalar field environment sourced by its companion. In this work, we calculate numerically the sensitivities of nonrotating black holes, including spontaneously scalarized ones, in three classes of Einstein-scalar-Gauss-Bonnet gravity: the shift-symmetric, dilatonic and Gaussian theories. When possible, we compare our results against perturbative analytical results, finding excellent agreement. Unlike their general relativistic counterparts, black holes in Einstein-scalar-Gauss-Bonnet gravity only exist in a restricted parameter space controlled by the theory's coupling constant. A preliminary study of the role played by the sensitivities in black hole binaries suggests that, in principle, black holes can be driven outside of their domain of existence during the inspiral, for binary parameters which we determine.

Quasinormal modes of slowly-rotating black holes in dynamical Chern-Simons gravity

The detection of gravitational waves from compact binary mergers by the LIGO/Virgo Collaboration has, for the first time, allowed us to test relativistic gravity in its strong, dynamical, and nonlinear regime, thus opening a new arena to confront general relativity (and modifications thereof) against observations. We consider a theory which modifies general relativity by introducing a scalar field coupled to a parity-violating curvature term known as dynamical Chern-Simons gravity. In this theory, spinning black holes are different from their general relativistic counterparts and can thus serve as probes to this theory. We study linear gravito-scalar perturbations of black holes in dynamical Chern-Simons gravity at leading order in spin and (i) obtain the perturbed field equations describing the evolution of the perturbed gravitational and scalar fields, (ii) numerically solve these equations by direct integration to calculate the quasinormal mode frequencies for the dominant and higher multipoles and tabulate them, (iii) find strong evidence that these rotating black holes are linearly stable, and (iv) present general fitting functions for different multipoles for gravitational and scalar quasinormal mode frequencies in terms of spin and Chern-Simons coupling parameter. Our results can be used to validate the ringdown of small-spin remnants of numerical relativity simulations of black hole binaries in dynamical Chern-Simons gravity and pave the way towards future tests of this theory with gravitational wave ringdown observations.

Quasinormal modes and late-time falloff of Finslerian black holes with cosmological constant

Quasinormal modes of scalar and electromagnetic field perturbations in Finslerian Reissner-Nordstr\"om black holes with a cosmological constant are investigated in this paper. We analyze the fundamental frequencies and dynamical evolution of quasinormal modes using the WKB approximation and finite difference method, respectively. Both approaches show that the periods of oscillation of quasinormal modes increase with higher Finslerian parameter ${\ensuremath{\epsilon}}^{2}$ if the multipole quantum number $l\ensuremath{\ge}2$. Using the Prony method, we show that the results obtained from the finite difference method are consistent with the results obtained from the WKB approximation. Quasinormal modes of Finslerian black holes possess spectrum splitting, which reflects the fact that the spherical symmetry of the Finslerian black holes is broken. The effects of the Finslerian parameter ${\ensuremath{\epsilon}}^{2}$ on the late-time tails of scalar and electromagnetic field perturbations are shown. The late-time tails of the Finslerian Reissner-Nordstr\"om black holes possess a power-law falloff. The power-law index has a discontinuous jump, while the Finslerian parameter ${\ensuremath{\epsilon}}^{2}$ varies from 0 to nonzero. Such a fact reflects the fact that asymptotic-infinity behaviors of the Finslerian Reissner-Nordstr\"om black holes are different from their counterparts in general relativity.

Effect of screening mechanisms on black hole binary inspiral waveforms

Scalar-tensor theories leaving significant modifications of gravity at cosmological scales rely on screening mechanisms to recover general relativity (GR) in high-density regions and pass stringent tests with astrophysical objects. Much focus has been placed on the signatures of such modifications of gravity on the propagation of gravitational waves (GWs) through cosmological distances while typically assuming their emission from fully screened regions with the wave generation strictly abiding by GR. Here, we closely analyze the impact of screening mechanisms on the inspiral GW waveforms from compact sources by employing a scaling method that enables a post-Newtonian (PN) expansion in screened regimes. Particularly, we derive the leading-order corrections to a fully screened emission to first PN order in the near zone, and we also compute the modifications in the unscreened radiation zone to second PN order. For a concrete example, we apply our results to a cubic Galileon model. The resulting GW amplitude from a binary black hole inspiral deviates from its GR counterpart at most by one part in ${10}^{2}$ for the modifications in the radiation zone and at most one part in ${10}^{11}$ due to next-order corrections to the fully screened near zone. We expect such modifications to be undetectable by the current generation of GW detectors, but the deviation is not so small as to remain undetectable in future experiments.

Compact stars in the Einstein dark energy model

We investigate the properties of high density compact objects in a vector type theory, inspired by Einstein's 1919 theory of elementary particles, in which Einstein assumed that elementary particles are held together by gravitational as well as electromagnetic type forces. From a modern perspective, Einstein's theory can be interpreted as a vector type model, with the gravitational action constructed as a linear combination of the Ricci scalar, of the trace of the matter energy-momentum tensor, and of a massive self-interacting vector type field. To obtain the properties of stellar models we consider the field equations for a static, spherically symmetric system, and we investigate numerically their solutions for different equations of state of quark and neutron matter, by assuming that the self-interaction potential of the vector field either vanishes or is quadratic in the vector field potential. We consider quark stars described by the MIT bag model equation of state and in the Color Flavor Locked phase, as well as compact stars consisting of a Bose-Einstein Condensate of neutron matter, with neutrons forming Cooper pairs. Constant density stars, representing a generalization of the Interior Schwarzschild solution of general relativity, are also analyzed. Also, we consider the Douchin-Haensel (SLy) equation of state. The numerical solutions are explicitly obtained in both standard general relativity, and the Einstein dark energy model and an in depth comparison between the astrophysical predictions of these two theories are performed. As a general conclusion of our study, we find that for all the considered equations of state a much larger variety of stellar structures can be obtained in the Einstein dark energy model, including classes of stars that are more massive than their general relativistic counterparts.